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Abstract. Recently, Imambekov and Glazman [Phys. Rev. Lett. 100, 206805 (2008)] showed 
that the dynamic structure factor (DSF) of the ID Rose gas demonstrates power-law behaviour 
along the limiting dispersion curve of the collective modes and calculated the corresponding 
exponents exactly. Combining these recent results with a previously obtained strong-coupling 
expansion we present an interpolation formula for the DSF of the ID Bose gas. The obtained 
expression is further consistent with exact low energy exponents from Luttinger liquid theory 
and shows nice agreement with recent numerical results. 



Cigar-shaped traps with cold alkali atoms have recently been used to obtain a quasi- ID 
quantum degenerate Bose gas, where atomic motion in the transverse dimensions is confined to 
zero-point quantum oscillations, in weak and strong interaction regimes [HIS]- Theoretically, we 
may describe the system as a one-dimensional rarefied gas where interactions of bosonic atoms 
can be described well by effective J-function interactions [3]. Thus the Lieb-Liniger model 
[H [5] is applicable. Being exactly solvable in the uniform case, the model, however, does not 
admit complete analytic solutions for the correlation functions. Up to now, this has been an 
outstanding problem in ID physics [6l[7]. Here, we propose an approximate formula for the DSF 
of the Lieb-Liniger gas that is consistent with known results in accessible limits and power laws. 

Dynamical density-density correlations, which can be measured by the two-photon Bragg 
scattering [81 [9], are described by the dynamic structure factor (DSF) [TO] 
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Here, we introduce the density fluctuations 5p{x, t) = p{x, t) — n and the equilibrium density of 
particles n = N/L. We consider the case of zero temperature, where (0| . . . |0) means ground- 
state average. The DSF is proportional to the probability of exciting the collective mode from 
the ground state with momentum q and energy hu> transfer, as one can see in the energy 
representation of Eq. ([1]) 

S{k,Lo) = Y.m5pk\n)\H{hu; - En + Eo), (2) 
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where 5pk = Y^j e~^^^i is the Fourier component of 5p{x). 



The Lieb-Liniger model [H [5] represents a uniform ID system of spinless bosons of mass 
m, interacting with pairwise point interactions V{x) = gBS{x); the interaction strength is 
assumed to be positive. Periodic boundary conditions are imposed on the wave functions. The 
strength of interactions can be measured in terms of the dimensionless Lieb-Liniger parameter 
7 = mg-Q/ {fi^n). Within the Lieb-Liniger model, the DSF has the following well-established 
properties. 

i) Luttinger liquid theory predicts a power-law behaviour of the DSF at low energies in 
the vicinity of the momenta = 0, 27rn, Aim . . . and yields model-independent values of the 
exponents [IH I12j. In particular, one can show |12l I13j that in the vicinity of "umklapp" point 
{k = 27rn, to = 0) 

S{k,uj) ^ {to^ -u;'i)^-\ (3) 

where K = hTrn/{mc) and c is sound velocity. Furthermore, within the Luttinger-liquid theory, 
the dispersion is linear in vicinity of the umklapp point: uj^[k) ~ c|/c — 27rn|. Relation ([3]) leads 
to different exponents precisely at the umklapp point and outside of it: 

afi N f k = 27rn, 



ii) By using in a non-trivial manner the Bose- Fermi mapping in ID [14], the authors developed 
the time-dependent Hartree-Fock scheme [151 [16] in the strong-coupling regime with the small 
parameter I/7. The scheme guarantees validity of the DSF expansion |15[ 116] 
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for u)- < 00 < a;_|_, and zero otherwise. Here uj±{k) are the limiting dispersions that bounds 
quasiparticle-quasihole excitations [5]. In the strong-coupling regime they take the form 
uj^(k) = h\2k-pk ± A;2|(l - 4/7)/(2m) + O (7-2). By definition, k-p = nn and ep = Tt^kl/{2m) 
are the Fermi wave vector and energy of a non-interacting Fermi gas, respectively. 
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Figure 1. Limiting dispersions uj± versus 
wave vector k for the coupling parameter 
7 = 5. The data are obtained numerically 
by solving Lieb-Liniger's system of integral 
equations [5l [6] . 
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Figure 2. Typical behaviour of the exact 
exponents in Eq. ([6|). The plot shows //± 
for 7 = 5 obtained numerically using the 
method of Ref. [IT]. 
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Figure 3. The Dynamic Structure Factor (DSF) in the thermodynamic limit. The proposed 
approximation ([8]) (hne) is compared to numerical data from Caux and Calabrese [19] (open 
dots). The dashed (red) line shows the data of Eq. ([8]) convoluted in frequency with a Gaussian 
of width Au;\/2 In 2 = 0.07eF/?i in order to simulate smearing that has been used in generating 
the numerical results of Ref. [19]. The numerical data of Ref. [19] suggests that contributions 
from multi-particle excitations for uj > uj^ (sharp line in parts a and b) are very small. Such 
contributions are not accounted for by the formula ([8]). 



iii) As was shown by Imambekov and Glazman [T7j, in the Lieb-Liniger model the DSF 
demonstrates power- law behaviour near the borders iOj-{k) 
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The positive exponents //j- [18] are related to the quasi-particle scattering phase and can be 
easily evaluated by solving a system of a few integral equations in thermodynamic limit |17| . 
We obtain the exact relation 
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which obviously differs from the Luttinger liquid exponent (j4]) for k ^ 27rn. However, 
Imambekov's and Glazman's result ([7]) is correct in the immediate vicinity of uj± provided that 
the finite curvature of uj-.{k) is taken into consideration. Thus the difference in the exponents 
can be treated [17] as an artifact of the linear spectrum approximation in the Luttinger liquid 
theory. Note, however, that the thin "strip" in uj-k plane where the exponents are different 
vanishes in the point k = 27rn; hence, the Luttinger exponent 2{K — 1) should be exact here. 

iv) The DSF can be calculated numerically by means of algebraic Bethe ansatz |19j . 

Here we suggest a phenomenological expression, which is consistent with all of the above- 
mentioned results. It reads 

Sik,-)=Ch^^-^ (8) 
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for LU-{k) < LV < u>-^-{k), and zero otherwise. Here C is a normalization constant, fJ'+{k) and 
H-{k) are the exponents of Eq. (I6|), and a = 1 + \/^/K. The normalization constant depends on 
momentum but not frequency and can be determined from the /-sum rule (see, e.g., Ref. |10j ) 
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We assume that in Eq. ([8|) the value of the exponent H-{k = 2im) comcides with its Umiting 
value d?]) in vicinity of the umklapp point. 
Now it can be easily seen from ([8]) that 



Thus, the suggested formula is consistent with the both the Luttinger liquid behaviour at the 
umklapp point and Imambekov's and Glazman's power-law behaviour in vicinity of it, as it 
should be. 

In the strong-coupling regime, Eq. ([8]) correctly yields the first order expansion Q. In order 
to prove this, it is sufficient to use the strong-coupling values of K = 1 -|- 4/7 -|- 0(1/7^), 
fi±{k) = 2 arctan[A;/(n7)]/7r -|- 0(1/7^) (see Ref. j20]), and the frequency dispersions. 

Comparison with numerical data by Caux and Calabrese |19j (figure E]) shows that the 
suggested formula nicely works in the regimes of both weak and strong coupling. 

Concluding, we propose the approximate formula dH) for DSF of the one-dimensional Bose 
gas at zero temperature. It neglects, in effect, only the contributions of multiparticle excitations 
outside the bounds given by the dispersion curves uj±, whose contribution is small. Our formula 
is consistent with predictions of Luttinger liquid theory, has the exact exponents at the edge of 
the spectrum, gives the correct first-order expansion in the strong-coupling regime, and shows 
nice agreement with available numerical data. 
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